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Aims. In this paper we compare the performance of multi and single-mode interferometry for the estimation of the 
phase of the complex visibility. 

Methods. We provide a theoretical description of the interferometric signal which enables to derive the phase error in 
presence of detector, photon and atmospheric noises, for both multi and single-mode cases. 

Results. We show that, despite the loss of flux occurring when injecting the light in the single-mode component (i.e. 
single-mode fibers, integrated optics), the spatial filtering properties of such single- mode devices often enable higher 
performance than multimode concepts. In the high flux regime speckle noise dominated, single-mode interferometry 
is always more efficient, and its performance is significantly better when the correction provided by adaptive optics 
becomes poor, by a factor of 2 and more when the Strehl ratio is lower than 10%. In low light level cases (detector noise 
regime), multimode interferometry reaches better performance, yet the gain never exceeds ~ 20%, which corresponds to 
the percentage of photon loss due to the injection in the guides. Besides, we demonstrate that single-mode interferometry 
is also more robust to the turbulence in both cases of fringe tracking and phase referencing, at the exception of narrow 
field of views (£> 1 ). 

Conclusions. Our conclusion is therefore that, from a theoretical point of view and contrarily to a widespread opinion, 
fringe trackers built using single-mode optics should be considered as a solution both practical and competitive. 

Key words. Instrumentation: interferometers - Methods: analytical - Techniques: interferometric 



1. Introduction 

Performance of ground-based optical interferometers is 
severely limited by the atmospheric turbulent piston which 
introduces a random optical path difference between the 
beams that are combined to produce fringes. These fringes 
are thus randomly moving on the detector, blurring the 
signal and preventing to integrate on time longer than the 
coherence time of the atmosphere, typically a few tens of 
milliseconds in the infrared. As a consequence, the limiting- 
magnitude and ultimate precision of ground-based interfer- 
ometers are dramatically reduced. 

By estimating and compensating in real-time the interfer- 
ometric phase - in other words by "locking" the fringes on 
the detector, fringe tracking and phase referencing devices 
are powerful instruments to circumvent this problem, al- 
lowing to integrate the interferograms on much longer time 
frames. Such instruments, which noticeably improve the 
sensitivity and the accuracy of interferometers are since a 
few years undergoing major developments, as can attest the 
amount of conce pts which are currently s t udied for various 
interferomet ers (|Le Bouquin et ah] |2008 
VLTI; |Berger et al.||2008| CHATTA; 
MROI). 



2008 



2008 



Sahlmann et al. 
Jurgcnson et al. 



In all the different concepts proposed, one key issue still 
in debate among the instrumental community is the rele- 
vance of using spatial filtering devices such as single-mode 
fibers or integrated optics to carry /combine the beams. It 
is often argued that spatial filtering is not suited for that 



type of instruments since only a fraction of the total flux 
is injected in the spatial filter component. This coupling 
efficiency, which is typically of the order of the Strehl ratio 
( Coude du Foresto et al.pOOO l, can indeed be low in pres- 
ence of strong turbulence and/or poor adaptive optics (AO) 
correction. However, the issue is not that simple. Single- 
mode filters only keeps the part of the incoming flux which 
is related to the coherent part of the corrugated wavefront. 
In other words single-mode devices only propagate the first 
mode of the electro-magnetic field leaving at its output a 
plane wavefront - that is a deterministic signal, the price 
to pay being a loss of flux correlated to the strength of 
the turbulence. At the contrary, multimode schemes such 
as bulk optics are preserving the total flux but keeps at the 
output the incoherent part of the wavefront (producing ran- 
domly moving speckle patterns in the image) which makes 
the interferogram sensitive to turbulence. Hence choosing 
between single-mode and multimode schemes corresponds 
to decide whether losing flux or losing coherence in the sig- 
nal will be the best strategy to optimize the performance 
of the interferometer. 

As a matter of fact, the efficiency of single- mode devices in 
terms of precision and robustness of the estimation of the 
amplitude of th e visibility has already been de monstrated 
experimentally (Coude du Foresto et al. 1997) and theo- 



retically as well ( Tatulli et al. |2*004 ) . However concerning 
the estimation of phase - which is the quantity of interest 
for fringe tracking and phase referencing instruments, the 
situation is not clear yet. If simulations have been initiated 
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Fig. 1. Sketch of the effect of the turbulence in interferom- 
etry. On each telescope pupil P(r) separated by the baseline 
distance Bi 2 , the wavefront is affected by a random phase 
screen, respectively 4>i(r) and (j> 2 (r). The strength of the 
turbulence depends on the parameters D/ro where D is 
the diameter of the telescope and ro is the Fried parameter 
(Fried 1965), describing the typical size of a coherent cell 
in the wavefront. Furthermore, the optical path difference 
between the two beams is randomly shifted by the turbu- 
lent piston <p\ 2 which is the phase difference between the 
two telescopes, averaged on the pupils. If an AO system is 
present on each telescope, then the wavefront is partially 
corrected in real time, and the quantity of interest are the 
residual phases (j)l(r) and (j) 2 (r) and the associated residual 
turbulent piston. 



in som e specific cases to es timate the ef fects of spati al fil- 
tering ( Buscher et al.|2 008: phase jumps; Tubbs|2 005 



com- 

putation of the coherence time), no theoretical formalism 
regarding this issue has, to our knowledge, been presented 
so far. 

In this paper, we derive the error of the interferometric 
phase in presence of detector, photon and atmospheric 
noises, both for the single-mode and multimode cases in 
presence of partial AO correction, following a formalism 
analogue t o the one previously developed for the squared 

p004| 



visibility in Tatulli et al. 



We then compare the per- 
formance of single-mode and multimode schemes for the 
estimation of the phase of the interferograms, and apply 
our analysis in the framework of fringe-tracking and phase 
referencing methods. 

2. The phase in ground-based interferometry 

2.1. General description and underlying assumptions 

The ground-based interferometer: As sketched in Fig. 
[T] we consider an interferometer with two telescopes, sep- 
arated by the baseline distance B i2 . We assume that both 
telescope are identical, that is they are described by the 
same pupil function P(r) of diameter D, therefore by the 
same collecting area £p defined as £p = J[P(r)] 2 dr. Their 
transmission are however different, namely t\ and t 2 re- 
spectively. The wavefront over each telescope is affected 



by randomly moving phase screens, respectively </>i(r) and 
(j) 2 (r). The typical cell size of these turbulent phases is ro 
( |Fried 1965 ) and the strength of the turbulence is character- 
ized by the quantity D/r$. When Adaptive Optic systems 
are present on each telescope, they partially compensate 
in real-time the wavefront corrugations of the atmosphere. 
The correction is however not perfect and it remains resid- 
ual phases 4>\(r) and (f> 2 (r) affecting the two telescopes. 
Because of this turbulence, the optical path difference (opd) 
between the two light paths to the star through the two tele- 
scopes is randomly shifted by the so-called turbulent piston 
</>i 2 , i.e. the phase difference between the two telescopes av- 
eraged over the pupils: 



%2 



J P(r ) 



r)dr — / P(r)(j) 2 (r)dr 



The interferogram: A 2-telescope interferogram consists 
in an incoherent part, which is the sum of the photomet- 
ric fluxes coming from both telescopes, and a modulated 
coherent part - namely the fringes, proportional to the so- 
called complex coherent flux F£ 2 which depends on the flux 
N (photons per surface unit and per time unit) and on the 
complex visibility Via of the source corresponding to the 
baseline frequency fi 2 — B\ 2 j\ of the interferometer. A is 
the effective wavelength of the interferogram. The modula- 
tion of the fringes can be whether temporal when the opd 
is scanned with moving piezo-mirrors (temporal coding), 
whether spatial when the opd is scanned with dedicated 
output pupils which separation defines the frequency of the 
modulation (spatial coding). 

Finally, the phase of the interferogram <& 12 , that is the 
phase shift with respect to the zero opd, originates from 
the source phase <fii 2 , a potential instrumental phase (j^yf , 
that we will assume equal to zero in the following] and the 
turbulence piston phase §\ 2 . 

2.2. Estimating the interferometric phase and its associated 
error 

Estimating the interferometric phase $12 requires to de- 
fine and compute from each interferogram an appropriate 
complex estimator F£ 2 of the complex coherent flux Ff 2 . 
Several methods are possible to build this estimator such 
as "ABCD" techniques in the image plane (Colavita"ll999 



Tatulli ct al. 2007) or analysis in the Fourier p lane (Roddier 
|fc Lena||1984| |Coude du Foresto et al. 1997). In any case, 
the measure' 
estimator F£ 2 : 



d phase is then the argument of the complex 



$12 = atan 



Im F- 



Re 



(2) 



In this framework, one can show that in first approxima- 
tion, i.e. when the error on F£ 2 is small compared to its am- 
plitude, the variance of the instantaneous phase (i.e. mea- 
sured over one interferogram) can be expressed as following 



1 considering a non zero instrumental phase (^{^ would only 
introduce an extra shift of the interferogram but would not affect 
the performance of the interferometer, as long as this instrumen- 
tal phase is stable or calibratable. 
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Table 1. Expression of the variance of the phase for the three noise regimes, respectively detector, photon and atmospheric 
regimes, both in the multimode and single-mode cases. See text and Appendix [X] for a description of the parameters 
involved. 
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( Chelli(l989l see Eq. (5), case j = k): 
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Re 
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(3) 



where E denotes the expected value. Chelli ( 1989 ) has also 
demonstrated that the variance of the interferometric phase 
is independent of the object phase. For sake of simplicity, we 
will therefore consider that the source of interest is centro- 
symmetric, namely 4>°^ = 0. 

Regardless of the method used to build the estimator from 
the interferogram, the statistics of the phase will depends 
on whether the interferometer is using a multimode or a 
single-mode design. In the following, we explore how the 
atmospheric spatial fluctuations of the turbulent wavefront 
affects the fringe pattern, in both multimode and single- 
mode cases. 



2.3. The coherent flux in multimode interferometry 

In multimode interferometry, the total flux on the detector 
remains constant (neglecting scintillation) and the interfer- 
ograms are created in speckle patterns randomly moving 
with the fluctuations of the turbulent wavefronts over the 
two telescopes. 

There are two different ways to combine the beam s in mul- 
timode interferometry (Chelli & Mariotti 1986): whether 
in the image plane, a technique which is suited to per- 
form spati al coding of the fringes on the detector, as known 
as Fizeau (Beckers & Hege 1984) and Michelson (Mourard 



et al.|l99~4 1 mountings^] or in the pupil plane which is com- 
monly used to temporally sample the inte rferogram throug h 
dedicated moving piezo-electric mirrors dPycket al.(l995] |. 
In both cases however, the expression of the coherent flux 
can be written as following (as shown in the demonstra- 
tion of Appendices |A.3.1| and |A.3.2 of this p aper and in 
Appendices Al and A2 of |Buscher et al.||2008"l ): 



F? 2 = S P iV A 



(4) 



T12 is the normalized interferometric transfer function re- 
sulting of the cross-correlation of the pupil of the two 
telescopes, corrugated by the residual atmospheric screen 
phases <j)\(r) and 02( r )j ( see Fig. [TJ. It writes: 



2ia 



J , [P( r )]2gi(05;(r-)-^(r)) fa 



(5) 



As an analogy wit h classical AO systems on single-pupil 
telescopes (see e.g. Fusco & Conan 20041, |Ti 2 | 2 can be 
seen as the instantaneous interferometric Strehl ratio, in 
the case of multimode interferometry. 

2.4. The coherent flux in single-mode interferometry 

The main property of single-mode devices such as fibers or 
integrated optic chips is to perform a spatial filtering of the 
input wavefront so that only its Gaussian part is transmit- 
ted. As a consequence, a single-mode device turns the input 
spatial wavefront fluctuations into intensity fluctuations at 
the output. Doing so, each outgoing wavefront is flat, hence 



2 in Fizeau combination, the output pupil is homothetic to the 
entrance one, whereas in Michelson scheme this property is not 
conserved. 
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Fig. 2. Global performance comparison between multimode and single- mode interferometry, all sources of noise being 
considered. Top: Error of the phase in multimode (dotted lines) and single-mode (solid lines) as the function of the 
number of detected photoevents (per time unit) and K-band magnitude. For both cases, we can see the three regimes: 
the detector noise (in 1/JV). then (shortly) the photon noise (in 1/ \/jV) and finally the speckle noise (saturation). Bottom: 
Ratio of multimode vs. single-mode phase error. Plots are shown for several levels of correction (Strehl ratios). The fixed 
parameters are adet = 4e~/pix, \Vij\ — 1, N p i X — 4. Two turbulence strengths are considered D/ro = 8 (left) and 
D/rg — 2 (right) which correspond to the average turbulence conditions of the VLTI with UTs (D = 8m) and ATs 
(D = 1.8m) respectively. For the conversion between detected photoevents and magnitude scale, we have considered a 
total transmission coefficient of 15%, a detector quantum efficiency of 50%, an integration time of 30ms and a spectral 
resolution of 35 with 2 pixels by spectral channel. 



the shape of the interferogram is deterministic, only de- 
pending on the instrumental configuration. The trade-off is 
that only a fraction of the flux is transmitted, correspond- 
ing to the coherent energy of the turbulent wavefront. The 
single-mode in stant aneo us coh erent flux thus takes the form 
(see Appendix pP} Eq. (|A.38|): 



^2 



(6) 



where p 12 (V) is the interferometric coupling coefficient that 
depends both o n the source extensi on and on the level of 
AO co rrection (| Tatull i et al.||2004[ see Eq. (3), and Eq. 



(A. 35 1, Appendix A.4| of this paper). Focusing on compact 



sources, that is astrophysical objects unresolved by a sin- 
gle telescope, p 12 (V) takes a simple express ion of the form 



p 12 (V) = p V 12 pi2 (Tatulli k Chelli 2005 Eqs. (4,5,6)) 
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The coherent flux thus rewrites: 



F L 12 = 



/ tlhpoVl2Pl2 



where p\2 is independent of the source properties: 
J p( r )e l WM)dr / P(r)e-^ r »dr 



P12 



J P(r)dr 



JP(r)dr 



(7) 



(8) 



and where po is the maxi mum achievable coupling effi- 
ciency, shown to be - 80% ( |Shaklan fe Roddier||1988| ), due 



to geometrical mismatch between the telescope Airy disk 
profile and the Gaussian profile of the propagated mode. 
Eq. ^ has to be compared with Eq. We can see that 
they are almost similar, but with one major difference: in 
Eq. d5l), the product of the phasors is integrated over the 
pupilwhereas in Eq. (JsJ , each phasor is first integrated 
over the pupil, then the product is performed. As noticed 



by Buscher et al. (2008), such difference is the mathemati- 
cal expressTon~^f the^roperty of spatial filtering. 
Again, \pi2\ 2 - which is equal to 1 in the case of perfect 
AO correction/absence of atmospheric turbulence - can be 
seen by as the instantaneous interferometric Strehl ratio, 
but for the single-mode case. 

3. Performance comparison 

3.1. Phase noise 

From Eqs. ([3jl4j7]), we derive the expression of the error 
of the phase for both single-mode and multimode cases, as 
developed in Appendix [5] We show that the variance of the 
interferometric phase can be decomposed as the quadratic 
sum of three terms, corresponding to the three regimes re- 
spectively detector, photon and atmospheric regimes: 



2 

°det d 



phot^ 



atm^ 



(9) 



where the detail of each term is given in Table [T] The phase 
error depends on parameters coming from the source and 
the instrument which we recall here: N is the number of 
incoming photons per surface unit and per time unit, | V12 1 
is the amplitude of the visibility, t\ and t% are the trans- 
missions of both telescopes, N P i x is the number of pixels 
to sample the intcrferogram, adet is the detector noise. The 
phase error also depends on atmospheric terms: the long 
exposure Strehl ratio S and its associated coherent energy 

e a *r where a , is the variance of the residual phase of 



a single telescope (Noll 1976), and of the phase structure 
function in presence or partial AO correction T)^ T {r) as 
modelled in Appendix |D) 

Figure [2] (top) shows the multimode and single- mode phase 
error as a function of the number of incoming photoevents 
(alternatively, source K-hanA magnitude), for different lev- 
els of AO correction (i.e. different Strehl ratios). In a general 
way, the behavior of the error is in both cases similar: at 
low fluxes, the detector noise dominates with a slope of the 
error in 1/N, then it goes by a (short) photon noise regime 
in 1/ v^/V and eventually reaches for the brightest sources 
a plateau due to speckle noise - analog to the one known 



in th e case of the visibility (Goodman 1985 Tatulli et al. 
20041 - which limits the ultimate precision on the phase. 
Note already that in order to keep the phase noise below 
reasonable levels, that is smaller than 1 radian, it is manda- 
tory to equip interferometers with AO systems, to at least 



insure a low-order correction of the turbulent phase provid- 
ing Strehl ratios greater than ~ 0.1 for -D/ r o = 8 and ~ 0.5 
for D/r = 2. 

Going in further details, Fig. [2] (bottom) shows the ratio of 
the error of the multimode phase by the error of the single- 
mode phase as a function of the number of photoevents. 
Clearly, two different behaviors occur whether we are con- 
sidering photon-starved or photon rich regime, as discussed 
in the following. 



3.2. Low light level regime - limiting magnitude 

In the case of faint sources, the phase error is dominated 
by the detector noise Udet^- The comparison between 
both multimode and single-mode cases is however not 
straightforward because the number of pixels required 
to sample the fringes is specific to each technique, and 
depends on the chosen combination scheme. Let us first 
review here the possible technical solutions: 
Multimode - (image plane) multiaxial combination: In this 
case, each speckle in which the interferograms are formed 
must be sample d correctly, that is it mus t be crossed by at 
least two pixels ( |Chelli fc Mariotti||1986 



the whole image of size A/r , the tota 



Then, if we take 
number of pixels 



2D/ro- As a consequence, in such 
N pix is dependent of the turbulence. 



required becomes N p i 
a combination mode, 

In the case of the VLTI at Cerro Paranal, the average ro is 
~ lm in the K-band, which gives turbulence strengths of 
D/r Q ~ 8 and D/r ~ 1.8 for the UTs (D = 8m) and ATs 
(D = 1.8m). This corresponds to a number of pixels that 
can vary quite a lot with the strength of the turbulence, 
respectively N pix = 16 (UTs) and N pix = 4 (ATs). 
Multimode - (pupil plane) coaxial combination: Providing 
that the interferogram is scanned faster than the coherent 
time of the atmosphere in order to "freeze" the fringes, the 
minimum number of pixels to obtain full information on 
the temporal interferogram is 3, which corresponds to the 
3 degrees of freedom: incoherent flux, fringe amplitude, 
and phase. However, instead of this "ABC" scheme with 
120 degrees between the three channels, it is frequently 
more practical to implement an "ABCD" scheme with 90 
degree phase shifts, hence requiring 4 pixels. Note that 
in the framework of phase tracking, it is also possible to 
only measure the sine component of the fringe with "AC" 
scheme with 180 degrees between the two channels. In any 
case, the number of pixels needed in a coaxial pupil plane 
combiner is between 2 and 4. 

Single-mode - multiaxial combination: On the contrary 
of the multimode/multiaxial combination, the number 
of pixels is here independent of the turbulence as the 
shape and frequency of the interferogram are fixed by the 
design of the beam combiner. Typically, the interferogram 
consists in a sinusoidal signal with a frequency defined by 
the separation of the output pupils (so-called the coding 
frequency), and where its amplitude is modula ted by the 
Gaussian envelop e of the single-mode device. | Tatulli fc"| 



LeBouquin (120061) have shown that the optimum number 
of pixels which respects the Shannon criterion (> 2 pixels 
per fringe) and prevents from an overlapping of the 
photometric and interferometric peaks in the Fourier plane 
is N plx = 10. 

Single-mode - coaxial combination: As in the coaxial 
multimode case, one just need here to sample the inter- 
ferogram with respect to the 3 degrees of liberty. As a 
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Fig. 3. Limiting flux as a function of the Strehl ratio (detector noise regime case), for turbulence strengths of D/ro = 8 
(left) and D/r = 2 (right). Plots are shown for an unresolved (V=l), fairly resolved (V=0.5) and fully resolved (V=0.1) 
source, respectively. Curves are going by pair, for single-mode (solid lines) and multimode (dotted lines) cases. The fixed 
parameters have the same values than previous figures. 



result, the number of pixels is again between 2 and 4. Note 
that instead of the usual temporal coding, fringes can be 
scanned simu ltaneously thanks t o "ABCD-like" integrated 
chip devices ( Benisty et al.|[2009 ). 

As a consequence, coaxial schemes - both in multimode 
and single-mode cases - appear more appropriate since 
they are using substantially less pixels than multiaxial 
ones. We remark however that these conclusions apply 
in the framework of fringe tracking where pair-wise 
combinations are favored. They may differ in the context 
of interferometric imaging instrument where an important 
number of baselines is involved. In the following we will 



10.0 



consider coaxial combiners with N p . 
to the standard "ABCD" sampling. 



4, corresponding 



As multimode and single-mode fringe trackers finally 
require the same number of pixels, a straightforward 
comparison of the expressions of the detector noise shows 
that multimode combiners will achieve in this regime 
slightly better performance by a factor of 1/po, because 
single-mode spatial filters can not transmit 100% of the 
flux for the geometrical reasons mentioned in previous 
sections. If we define the limiting magnitude]^] such as the 
error of the phase is equal to 1 rad, which - apart from 
very low AO correction levels - occurs in this detector 
noise regime, the corresponding limiting flux is given by: 
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Fig. 4. Phase speckle noise as a function of the Strehl ra- 
tio for both single-mode (solid line) and multimode (dotted 
line) cases, with D/r n = 8. The ratio between the multi- 
mode and the single-mode phase errors is plotted in dashed 
line (y-axis on the right side). 



K, 



multi, single] 



V 2 y ^pix^det 



(10) 



3 note that the definition of the limiting magnitude depends 
on the estimator chosen to measure the phase. In this paper, 
the interferometric phase is merely estimated at the baseline 
frequency / 12, th at is at t he top of the interferometric peak (see 
Appendices A. 3 and A.4I. If a different estimator is used, like 
integrating the nigh frequency peak over the frequency range 
\.fi2 - D/X , fi2 + D/X] (see e.g. |Roddier fc Lena|1984| |Mourard 



et al.|[l994 in the case of the squared visibility estimators) the 
expression, and therefore the value of the limiting magnitude 
will change accordingly. 



with the factor 1 or po depending whether we are using 
multimode or single-mode interferometers. As illustrated 
by Fig. |3j the gain in limiting magnitude for multimode 
combiners is ~ 0.25 magnitude. 



3.3. High light level regime: the speckle noise 

In the photon-rich regime, the speckle noise dominates and 
the error on the interferometric phase reaches a plateau, 
that is, is not dependent on the flux of the source anymore. 
As illustrated on Fig. [4j single-mode interferometers always 
provides in such a case lower phase error than multimode 
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ones, emphasizing the remarkable properties of spatial fil- 
tering of the turbulent wavefront by single-mode devices 
Such behavior has been already noticed by |Tatulli et al 
(2004) for the estimation of the squared visibility, show- 
ing that the so-called modal speckle noise of the visibility 
was, for a given AO correction, always smaller than classi- 
cal speckle noise of multimode interferometers. The concept 
of modal speckle noise can also be applied for the single- 
mode interferometric phase. However at the difference of 
the squared visibility for which the modal speckle noise is 
for a point source, the phase modal speckle noise always 
exists, independently of the size of the source. 
Note that the gain of using single-mode schemes is all the 
more important than the the level of correction is low. If 
for fairly good correction with Strehl ratio above 20%, the 
difference remains marginal with a factor ~ 1 — 1.5 between 
the single-mode and multimode phase error, the situations 
where bright sources are observed with low/none AO cor- 
rection will highly benefit of single-mode interferometers. 
In such cases the precision of the phase can increase by 
at least a factor of 2 and much more when using spatial 
filtering of the corrugated wavefront with a typical Strehl 
ratio below 10%. This is a typical counter-intuitive exam- 
ple where it is more profitable to loose a substantial part of 
the flux and keep only the coherent part of the perturbed 
wavefronts than conserving the whole flux at the price of 
introducing additional atmospheric noise. 

4. Application to fringe tracking 

In this section, we apply the formalism developed previ- 
ously in the context of on-axis and off-axis fringe tracking, 
that is when the phase is estimated and compensated in 
real time to stabilize the fringes on the detector. We focus 
on the relative performance of fringe tracking systems using 
whether multimode or single-mode schemes. 



4.1. Coherent integration 

For each interferogram, fringes are shifted from the zero op- 
tical path difference by the turbulent piston. This piston is 
corrected in real time by the fringe tracker which estimates 
the phase of the fringes and compensates the optical path 
difference in real time by moving dedicated mirrors. 
Centering the fringes in real time allows to perform coher- 
ent integration of the signal, that is to integrate on time 
scales much longer than the coherence time of the atmo- 
sphere. However, the estimated piston used for the opd 
correction is affected by a random measurement error e(t). 
As a result, the interferogram is not perfectly centered and 
is still moving with an excursion depending on the statis- 
tics of the noise. If the interferograms are integrated over a 
time long enough for the realizations of the random error 
e(t) to span all the range of its probability law, then, by 
rules of ergodicity, the visibility will be affected by a loss of 
contrasiQ which writes: 



1.0 



< v > t =< v > £ = Ve- 



in) 



4 For sake of simplicity, we assume here that the phase is com- 
pensated instantaneously, hence not taking into account the de- 
lay of the fringe tracking loop between the measurement of the 
phase and its correction. The problem of time delay, which is in- 
dependent of the multimode o r single-mode natur e of the fringe- 
tracking system, is treated in Conan et al. (2000 1. 
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Fig. 5. Long exposure visibility attenuation due to the 
atmospheric noise of the measurement of the phase, as a 
function of the Strehl ratio for multimode (dotted line) and 
single- mode (solid line) systems, with D/tq = 8. The ratio 
of the attenuation between single-mode and multimode case 
is plotted in dashed line (y-axis on the right side). 



where <> e is the ensemble average over the realizations 
of the noise, and is the variance of this noise as com- 
puted in Sect. |3.l[ Hence, for long integration times, the fact 
that the instantaneous turbulent phase is corrected only to 
the precision of its estimation, introduces an attenuation of 
the coherent flux (equivalently, of the visibility) by a factor 

2 

_ 

e t~ . Such bias therefore depends on whether single-mode 
or multimode fringe tracking is used. 

For high-light level regime in which fringe tracking is mostly 
expected to work, Fig. [5] shows the attenuation of the visi- 
bility as a function of the Strehl ratio for both single-mode 
and multimode systems. We can see that acceptable loss 
of contrast, typically <; 0.8, is achieved as soon as mod- 
erate AO corrections with Strehl 0.1 (for D/ro = 8) is 
provided. On the contrary, for lower performance of the 
AO system, the attenuation coefficient drops rapidly and 
the advantage of tracking the fringes is lost due to the un- 
compensated turbulent fluctuations of the phase over the 
pupils. 

One important conclusion to draw is that, for a given AO 
correction, the visibility attenuation induced by the noise of 
the phase is always smaller when using single-mode fringe 
tracking systems instead of multimode ones. In other words, 
the maximum achievable atmospheric contrast is always 
higher using single-mode fringe tracker than multimode sys- 
tems. This is especially true for low AO correction cases 
where the attenuation can becomes twice larger and more, 
as shown in Fig. [5j 

4.2. Phase jumps 

It is usually implicitly assumed that the phase estimated 
in interferometry is effectively the shift of the fringes with 
respect to the zero opd, namely the piston phase. By defini- 
tion, the piston phase shift is the average differential phase 
across the apertures, as illustrated on Fig. [I] and described 
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Fig. 6. Density probability of the multimode instantaneous Strehl ratio (/ = | T\ 2 1 2 , red lines) and single-mode one 
(/ = I/O12I 2 , black lines) for various levels of AO correction and turbulence strengths. Top: constant turbulence strength 
D/rO — 8 and the number of corrected Zernike are, from left to right, n z = 10, 15,30. Bottom: constant AO correction 
n z = 3 (tip-tilt only), with various turbulence strengths, from left to right D/ro = 8,4, 2. 



in Eq. ([I]) that can be rewritten: 
# 2 (*)= / Wxir,t) 



i(r,*)]dE 



(12) 



following the formalism of |Buscher et al. (2008 see Eq. (6)), 
where dS is an elemental within the aperture S. 
As first brought up by this author, the practical estima- 
tion 0x2(0 °f such piston phase is performed by taking the 
argument of the complex phasor averaged over the pupils, 
that is: 



<^ 2 (i) = axg( I e^^-^^dE 



in the multimode case, and: 



cj)P 2 (t) = arg ( / e^^dE / e^^dE 



(13) 



(14) 



in the single-mode one. 



By comparing Eq. ( 12 ) to Eqs. (13 14 ) we can see that the 

true piston phase §\ 2 and the estimated one 4>\ 2 are liter- 
ally not the same, unless when there is no phase fluctuation 
(i.e. no turbulence or a perfect correction of the wavefronts) 
across the apertures. 

When considering a small amount of fluctuations S<fii2(r, t) 
such as (j>\(r,t) — <f>2{r,t) = S(f>i2(r, t), a Taylor expansion 



y 12 



and 



h 1 ' 

£L2- 



will then de- 



2008) 



of Eqs. ( 13|14 ) can show that 
viate roughly as J s [(50i2(r, i)] 3 d£ (Buscher et al 
But if these fluctuations are large (typically above 1 ra- 
dian), there can be strong differences between the piston 
phase and the argument of the complex phasor. More pre- 
cisely, strong divergences are occurring when the complex 
phasor has a very small amplitude and eventually crosses 
the origin of the complex plane. In such cases, rapid phase 
jumps of the argument of the phasor can be experienced 
whereas these jumps are not seen in the piston phase. As a 
consequence the correction performed by the fringe tracker 
may be highly wrong and fringes potentially lost, especially 
if the science camera is workin g at a different wave length 
than that of the fringe tracker (Buscher et al. 2008[ ), as it 



is the case e.g. for the FINITO instrument (fring e trac king 



Le 



in H-band and correc tion in K-band, Gai et al. (2004) 
|Bouquin et~aI1 ( |2008[ )). 
it 5 therefore important to know the probability of 
phase jumps to occur in order to ensure an observa- 
tional/instrumental context in which these events are 
avoided as much as possibl e. So far, the conclusions about 
this point are not clear : if Tubbs ( 2005 ) implied that spa- 
tial fi l tering could be the source of such anomalies, |Buscher 



et al. (2008) on the contrary argues that AO correction and 
most of all spatial filtering help to reduce these singulari- 
ties. All the previous analysis were however based only on 
simulations. We propose here a theoretical analysis of this 
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phenomenon. We emphasize that our analysis focuses on 
the rate of occurrence of phase jumps, as defined above. 
We do not study the consequences of these phase jumps in 
terms of effective loss of fringe tracking. Such a causality 
will depend on the system (e.g. single wavelength vs. multi- 
ple spectral channels methods) used to practically estimate 
and correct the opd. Considering that a phase jump would 
systematically induce a failure in the fringe tracking system 
therefore corresponds to the worst case scenario. 
The fact that the coherent flux drops to zero because of 
turbulent phase fluctuation depends on whether the inter- 
ferometric transfer function |Tia| for multimode systems or 
the interferometric coupling coefficient \p\2\ in the single- 
mode case, drops to zero. We thus want here to establish 
the probability density of these quanti ties, and study how 
likely they can take very small values. |Canales &: Cagigal] 
( 1999 ) have studied the distribution of speckle statistics in 



presence of partial AO correction for mono-pupil telescopes. 
They have shown that the density probability dp{I) of the 
intensity I at the center of the image - that is by definition, 
the instantaneous Strehl ratio - follows a Rician statistics 
of the form: 



dp(I) 



1 

2a 2 



I + a< 
2a 2 



Jo - 



(15) 



where a 2 and a 2 depends on the first and second order 
moments of the real and imaginary part of the complex 
phasor describing the AO corrected turbulent wavefront. 
Note that such analytical definition of the density proba- 
bility implicitly assumes a time-constant ro to characterize 
the turbulence. This hypothesis is the limitation of our 
model as ro may actually vary on time scales of minutes 
and shorter, driving to brief episodes of very small ro that 
would trip up AO systems and lead to very small Strehl 
ratio for a brief period of time. However, modelling this 
effect requires more complex and heavy simulations of 
partially AO-corrected turbulence which are out out the 
scope of this paper. 

We have emphasized in Sect.[2]that |Ti 2 | 2 and |pi2| 2 repre- 
sent the instantaneous interferometric Strehl ratio, for the 
multimode and single-mode cases, respectively. Therefore, 
by straightforward analogy with an interferometric pupil 
instead of a monolithic one, we can show that / = | T\ 2 1 2 
and / = I/O12 1 2 are also following a Rician distribution. 
The parameters a 2 and a 2 relative to both instantaneous 
interferometric Strehl ratio can be directly derived from 
the formalism of previous section and their expressions are 
given in Table |B.1| of Appendix [B| 

Figure [6] shows the density probability of |T 12 | 2 and | /o 12 1 2 
for various levels of AO correction and turbulence strength. 

High-order AO correction: First, we can see that, 
as soon as moderate AO correction is applied - roughly a 
few tens of modes or a Strehl ratio <; 0.3, the probability 
density displays a bell shape, all the more narrow than the 
AO correction is high, and the probability that |Ti2 1 or 
| pi2 1 goes to very low values is null (typically, for n z = 15 
and D/rO — 8, the probability for the interferometric 
instantaneous Strehl ratio to be < 0.01 is < 10~ 9 , and 
< 0.05 is < 0.02%), emphasizing again the relevance of as- 
sociating high-performance AO systems to fringe tracking 
devices, the wavefront being spatially filtered or not. We 
note that the intensity mean value (< / >= 2er 2 + a 2 ) is 



systematically slightly higher in the single-mode case and 
the dispersion (erf = 2er 2 + 4a 2 a 2 ) slightly smaller. The 
difference is however not critical. 

Low-order AO correction: For low-order AO such 
as tip-tilt correction only, the situation is fairly different. 
The shape of the distribution has changed, peaking at 
zero in the case of strong turbulence with high values of 
D/ro 5. The probability to have a very small value 
of the intensity is thus significant, and events such as 
phase jumps are likely to happen. In other words, tip-tilt 
correction is a too low order correction to prevent from 
such events. 

Figure [7] displays the probability P(I) for the instanta- 
neous interferometric Strehl ratio to reach values close 
to zero (i.e. P(I < e) with e very small) for different 
strengths of the atmospheric turbulence. We can see that, 
for strong turbulence (D/ro *S 5), the probability P(I) < e 
is always higher in the single-mode case. In such cases, 
we expect phase jumps to occur more frequently with 
single-mode fringe tracker, as the likelihood that two 
speckles are simultaneously entering the single-mode fibers 
of the two telescopes is weak. However, this does not mean 
that phase jumps should not occur often with multimode 
fringe tracker too. As an example the probability to have 
instantaneous interferometric Strehl lower than 0.01 in 
the multimode case is ~ 78% for D/ro — 8, whereas it is 
100% in the single-mode one. In other words, regardless of 
the multimode versus single-mode issue, one should never 
consider to perform fringe tracking with big telescope 
apertures associated with the sole tip-tilt correction. 
For weak turbulence (roughly D/ro & 5, Fig. |6j bottom 
middle and right), the distribution starts to look Tike a bell 
shape again, with however non null probability to have a 
zero intensity. But as shown in Fig.[7]this probability is now 
higher in the multimode case, as only a few and big speck- 
les are present in the images. As a matter of fact, this case 
is the one treated by Buscher et al. ( 2008[ see Fig. 4, with 
D/ro = 4, tip-tilt correction) in his simulations, and we 
come to the same conclusion than his: in these conditions, 
spatial filtering enables to decrease the number of phase 
jump events. We however emphasize here that this state- 
ment is true only for cases of moderate turbulence strength. 

Application to the VLTI: We recall that at Paranal the 
average r is ~ lm in the K-band, which gives turbulence 
strength of D/ro ~ 8 and D/r ~ 1.8 for the UTs and ATs, 
respectively. UTs come with high-order Adaptive Optics 
systems (MACAO, |Arsenault et al.|2004 ) therefore one can 
equally choose multimode or single-mode fringe tracking 
schemes since (i) phase jumps are unlikely to occur in any 
case, and (ii) multimode and single-mode fringe trackers 
will provide the same robustness to eventual phase jumps. 
In the case of ATs where tip-tilt correction is currently 
provided, the likelihood to endure phase jumps is also 
very low. However to maximize the stability of their fringe 
tracking system, it still seems appropriate to accompany 
the ATs with single-mode spatial filtering devices as the 
probability to undergo such phase jumps remains roughly 
10 times lower (see Fig. [7j right). Alternatively, providing 
higher order AO correction to ATs will fix the phase jump 
issue. 
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Fig. 7. Left: probability for the instantaneous interferometric Strehl ratio I to be lower than the value e, in the case of 
tip-tilt correction only, for different strengths of turbulence D/r = 8,4,2 (respectively in blue, magenta and red), for 
the multimode (dotted line) and single-mode (solid line) cases. Right: Multimode versus single-mode probability ratio. 




Aoc [arcsec] Aa [arcsec] 

Fig. 8. Left: astrometric phase error due to atmospheric turbulence, as the function of the angular distance Aa between 
the source and the reference star, for both single-mode (solid lines) and multimode (dotted lines). Right: ratio of the 
multimode versus the single-mode astrometric error. Results are shown for various AO correction levels as described in 
previous Figures. We have assumed one single turbulent layer located at h = 10km. 



4.3. Phase referencing - astrometry 

The interest of phase referencing inst ruments such as 
PRIM A recently installed on the VLTI ( jPelplancke et al. 
2000 1 , is (i) to provide sub-microarcsecond precision as- 
trometry, allowing e.g. to detect the presence of a faint 
companion (extrasolar planet) around the central star 
(Launhardt et al. 2008), and (ii) to drastically increase 
the limiting magnitude of the interferometer by locking the 
fringes of the (possibly faint) science object on a simulta- 
neously observed bri ght off- axis reference st ar whose phase 
is used as reference (Sahlmann et al. 20081. Such method 



requires that the nearby reference star is close enough 



typically in the isoplanatic patch - in order to assume that 
both the wavefront of the source and the reference are iden- 
tically perturbed by the atmospheric turbulence. Strictly 
speaking, this assumption is not true as the optical paths 
of the two stars arc different, crossing different part of the 
atmosphere in the turbulent layers. As a result, the loss of 
correlation between the two wavefronts will drive to an at- 
mospheric noise on the astrometric phase which will lower 
the ultimate performance of this method. 
We call respectively <f) s , and 4>^ ah , the phase of the astro- 
physical target and the phase estimated from the off-axis 
reference source, located at an angular distance Aa of the 
science object, and where we assume for the sake of sim- 
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plicity a single turbulent layer located at an height h from 
the ground. The astrometric phase is then simply defined 
by: 

A$ ah = </>, 



lAaii 



and its associated error writes: 



2ai-2cov(^^ atl ) 



(16) 



(17) 



Equation 17 tells us that the noise of the astrometric 



phase increases as long as the correlation between the two 
turbulent wavefronts cov(4> s , 4>^ ah ) decreases, that is as 
long as the spatial distance Aah between both turbulent 
phase screens is getting higher. The quantitative effect 
of the separation between the star and its reference on 
the interferometric narro w-angle astrometri c error has 
been already studied by 



Shao & Colavital (fl992| and 



we refer the readers to tneir paper for further details 
In our analysis, we focus on the relative performance of 
astrometry between multimode and single-mode scheme, 
studying how spatial filtering of the wavefronts will impact 
on the astrometric error, as developed in Appendix |Cj 
For sake of simplicity we have assumed in the following 
one single turbulent layer located at ft. = 10km as it is 



the stronges t layer at the Cerro Paranal site (Masciadri 
et al. 1999). There is also a strong layer very close to 



the ground at around 20m (Martin et al. 2000) but that 
only marginally contributes to the phase decorrelation 
since the linear distance at stake is 500 times smaller 
than the one relative to the high altitude layer. Figure [8] 
(left) shows the evolution of the astrometric error as a 
function of the angular distance between the source and 
the reference, for both single-mode and multimode schemes. 

Large field of view: When the reference source is 
further than roughly ~ 1 from the astrophysical object of 
interest, we can see that spatial filtering allows to increase 
the precision of the astrometric phase. The improvement 
is better as the AO correction is low, reaching a factor 
of ~ 4 without AO correction, as emphasized by Fig. 

t (right). This argues in favor of using spatial-filtering 
ements in the design of astrometric instruments with 
large field of view (FOV) as it has been chosen for PRIMA 
(FOV ~ 30 ). Looking in further details, one can notice 
that the phase tends asymptotically towards a plateau, 
from an angular separation of ~ 10 which roughly corre- 
sponds to the isoplanatic angle of Paranal in the K-band2 
Indeed from this angle, wavefronts can be considered as 
uncorrelated and the astrometric error converges to v^cr^, 
which does not depend on the separation Aa any longer. 
As a consequence in this regime (i.e. Aa > 10 ), the 
astrometric phase error ratio can be approximated by the 
speckle phase error ratio, which has already been discussed 
in Sect. fO 



Narrow field of view: At the contrary, when wavefronts 
are still strongly correlated (Aa ^ 1 ), the performance 
of single-mode astrometry is slightly lower (by a factor 
~ 1 to ~ 1.5) than that of multimode astrometry, the 



Martin et al. 



(2000) which 



considering the measurements of 

have found an average isoplanatic angle of ~ 1.9 in th e visible, 
and recalling t hat it evolves with wavelength as A 1,2 (Shao & 
Colavita| 1992 1. 



error of the latter decreasing faster as the separation 
between the astrophysical source and the reference shrinks 
to zero. By smoothing the turbulent wavefronts across the 
apertures, spatial filtering is indeed lowering the effect 
of the strong correlation between the wavefronts. Such 
situation concerns the cases for which the reference is very 
close to the star, like GRAVITY the second generation of 
astrometric instruments of the VLTI ( |Gillessen et al.| |2006 , 
FOV ~2"). 

5. Summary and conclusions 

In this paper, we have provided a theoretical formalism 
that allows to derive the error of the interferometric phase, 
both in cases of single-mode and multimode intcrfcromctry. 
From these derivations, we have demonstrated that: 

— Contrarily to a widespread idea, losing flux by injecting 
the light into single-mode spatial filters is not a per- 
formance killer for estimating the phase. As the mat- 
ter of fact, single- mode interferometry provides better 
performance than that of multimode interferometry, un- 
less when the interferometer is working in the detector 
noise regime (faint sources). In such cases, multimode 
interferometry is slightly better, providing a phase error 
smaller by a factor po — 0-8 which is the maximum frac- 
tion of flux that can be injected in single-mode devices. 

— In cases of bright sources observations, spatial filtering 
is shown to be very efficient, especially when the AO cor- 
rection of the turbulent wavefronts is poor or absent. In 
such situations, the precision of the single-mode inter- 
ferometric phase is better than that of the multimode 
one by a factor of 2 and more when the Strehl is below 
10%. 

— Single-mode interferometry also proves to be more ro- 
bust to the turbulence to both locking and coherently 
integrating the fringes and providing a better astromet- 
ric precision when using phase referencing techniques, 
except for narrow field of views (FOV 1 ) . 

In conclusion, from a theoretical point of view and con- 
trarily to a widespread opinion, single-mode fringe tracking 
should be seriously considered as an advantageous technical 
solution. Furthermore, the astronomers should realize that 
the many gains of single-mode interferometry (flexibility of 
the solutions, robustness to the alignment, less optical ele- 
ments,...) may compensate significantly a modest and lim- 
ited loss of performance in the faint sources case compared 
to multimode interferometry solutions. This is all the more 
true that the detectors should soon evolve to reach photon 
counting capability. 
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Appendix A: Computation of the interferometric phase error 

A.l. General formalism 

It follows Goodman's approach ( Goodman|198 5), based on a continuous model of detection process, where the recorded 
signal (i.e. the interferogram) can be represented as: 



A 



d(%, y) = y~] H x -xk,y- yk) 



k=i 



and its Fourier transform: 



A" 



-2i-n(f x x k +f v y k ) 



(A.l) 



(A.2) 



k=l 



where the position (x n ) and the number of photoevents (per time unit) (K) are statistical processes with probability laws 
depending on the intensity distribution I( x). To take also into account the detector noise, one needs to add Gaussian 
additive noise (e) to the previous equation ( |Tatulli et al.||2004| [TaFulli fc Chelli|2005[ ): 



A' 



Qif) = D(f) + e(/) = £ e-^'^+f'^ + e(f) 



(A.3) 



fc=i 



Note that this formalism is also valid for a temporal coding of the interferogram, that is replacing the spatial position 

(xkiUk) by the (1-D) temporal sampling t, and similarly the spatial frequency / by the temporal one v. 

As Q represents the spectrum of the interferogram, the phase of the spectrum is merely the argument of this estimator: 



atan(Im(Q),Re((3)) 



(A.4) 



In this framework, Chelli (1989) has shown that in first approximation (i.e. small noise error), the variance of the phase 
can be expressed as following: 

2 lE(|Q| 2 )-Re[E(Q 2 ) 
^=2~ 



[E(Q)P 

then using Goodman's formalism described above, one can show that ( |Chelli|p89l |Tatulli et aT^2004| ): 

E(|Q| 2 ) = <K 2 \ l (f)\ 2 + K + N plx al t 

E(Q 2 ) = <K 2 i 2 (f) >$ 
E(Q) = < Ki(f) >* 



(A.5) 



(A.6) 
(A.7) 
(A.8) 



where K is the average number of photoevents per time unit, adet is the detector noise and N p i X is the number of pixel 
to sample the signal. i(f) is the normalized spectrum of the interferogram (i.e. i(0) — 1). 

<>$ denotes the expected value with respect to the atmosphere. How the previous expressions are unfolding depends 
on whether we deal with multimode or single-mode interferometry. 



A.2. Useful definitions 

In this section we introduce the concepts and notations that we will use in our formalism to derive the expression of the 
noise of the interferometric phase. Let us consider (see also Fig. [I] in the paper): 

— an interferometer made of 2 telescopes described by their pupil functions Pi{r), Pi{ r ) and their associated transmission 
y/ti, \ft2- Note that with such a definition P\ and Pi are respectively centered at the position r\ = A/i and T2 = A/2, 
and we introduce the quantity /12 = f% — fi as the baseline frequency of the interferometer. 

— (f>i(r) and <t>2{ r ) the partially AO-corrected atmospheric turbulent phase screens. We also introduce T> ! p r the structure 
function of the residual phases (assuming the same level of AO correction for each telescope) together with ct^ , its 
associated residual phase variance (see also Appendix |d|) . 

— the incoming wavefront of the source Vt'(r). By definition we have ^(O)! 2 — N, N being the number of photons per 
surface unit and per time unit emitted by the source, and [^'(r) l I , *(r + A/)]/|5'(0)| 2 = V(f) the visibility of the 
source at the spatial frequency /. 

We emphasized that the following derivations will be done assuming monochromatic interferograms at the wavelength A. 
Indeed, taking into account a non-null spectral bandwidth 5X is equivalent to consider a monochromatic interferogram 
at the effective wavelength Ao of the filter, modulated in amplitude by an envelope, which width is fixed by the spectral 
coherence length L c = \q/5A. Moreover, such effect is independent of the multimode or single- mode nature of the 
interferogram. 
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A. 3. Phase noise in multimode interferometry 

In the following, we derive the estimator of the coherent flux of the interferogram, and subsequently the error of the inter- 
ferometric phase in the multimode case. We provide two formalisms that respectively describe both beam recombination 
schemes: in the image plane and in the pupil plane. 



= h 
+ h 



A. 3.1. Image plane recombination 

Image plane recombination is well suited to describe multiaxial interferometric schemes where the fringe pattern is 
spatially coded on the detector, such as for the GI2T and VEGA/CHARA instruments. Note that co-axial (i.e. temporal) 
coding can also be performed in the image plane but in practice - to our knowledge - no interferometric instruments are 
using this technique. 

In such case, the spatial distribution of the complex amplitude of the radiation field E(r) can be written as following: 

E{r) = VhPi(r) exp 4 ^ #(r) + >/^P 2 (r) exp^ (r) #(r) (A.9) 

By rules of diffraction theory, the Fourier Transform of the interference pattern formed in the image plane is the auto- 
correlation of the complex amplitude E(r), namely: 

/(/) = J E(r)E*(r + Xf)dr 

J Pi(r)Pi(r + A/) exp'MW-Ktr+W) *( r )**(r + Xf)dr 

i f P2{r)P 2 {r + A/) exp l «^ r >-^(*-+\f)) #( r )#*( r + A/)dr 

+ Vhh J Pi(r)P 2 (r + Xf)exp i ^ r) -^ r+Xf ^^(r)^f*(r + Xf)dr 

+ Vhh J Pi(r + A/)P 2 (r > )exp i( ^M-^('-+ A /))$*(r + A/)*(r)dr (A.10) 

Introducing |^(0) | 2 — N in previous equation enables to make the complex visibility function appear: 

/(/) = Nt 1 V(f)S 1 (f) + Nt 2 V(f)S 2 (f) + NVhhV(f)S 12 (f) + NVh^V(f)S* 12 (-f) (A.ll) 

with Si(f) and S 2 (f) the so-called photometric peaks resulting on the autocorrelation of each corrugated pupil, and 
812(f) the interferometric peak arising from the cross-correlation between both pupils: 

Si(f) = /'Pi(r)P i (r + A/)exp < WW-^^+A/)) dr) (i=[1)2 ]) (A.12) 

Si 2 (f) = J Pi(r)P 2 (r + A/) exp l WM-^( r+A/ » dr (A.13) 

From equation |A.11[ we can straightforward derive the number of photoevents: 

K = 1(0) = Z P N(t 1 +t 2 ) (A.14) 
which is turbulence independent (neglecting scintillation). Ep is the collecting area of a single telescope: 

E P = J [Pi(r)] 2 dr = J [P 2 (r)] 2 dr (A.15) 

assu ming for sake of simplicity that both pupils are identical. 

Eq. (A.ll) also shows that multimode interferometry with image plane recombination continuously transmits the whole 
spatial frequencies /. The complex visibility V12 — V(fi 2 ) is then derived from the estimation of the coherent flux at 
the baseline frequency of the interferometer fi 2 , that is: 

F? 2 = NVht2V 12 S 12 (f 12 ) (A.16) 

and the normalized spectrum of the interferogram takes the form: 

where Ti 2 (f) — S\ 2 (f) /T,p is the normalized interferometric transfer function, such as |r 12 (^12) I = 1. As a consequence, 
by analogy with the definition of the Strehl ratio for a single telescope, \Ti 2 (fi 2 )\ 2 can be seen as the instantaneous 
multimode interferometric Strehl ratio. 



E. Tatulli et al.: Multimode vs. single-mode interferometric phase, Online Material p 3 
A. 3. 2. Pupil plane recombination 



In such a schem e, fringes are formed in the pupil plane by means of geometrical operations of the entrance pupil ( Chelli & 
Mariotti||1986 1. This formalism is usually well suited to describe co-axial recombination with temporal coding where the 



pupils are superimposed at the beam splitter, hence the fringes formed in the pupil plane, as in the IOTA instrument. In 
practice, the fringes are form by translating one pupil ove r the other (t ranslation of vector /12), and by introducing an 
optical path delay (9i(t) and 9 2 {t)) on each beam (see e.g. Buscher||1988 Appendix A2) to modulate the fringe pattern. 
In this case, the complex amplitude of the superimposed electric fields writes: 

E(r) = VMMrOexp^i^ #(r) exp i9l(t) +VhP 2 {r + \f 12 ) exp l ^ (r+xfl2 *> *(r + A/ 12 ) exp ie2(t) 



(A.18) 

The fringe pattern I(t) is then formed by focusing the light on one single pixel of the detector. Again, by rules of image 
formation, it writes as following: 

(A. 19) 



I(t)= I \E(r)\Mr 
which, by making use of the expression of E(r) rewrites: 



1(f) 



Ej>JV(ti + i a ) + AVii^exp lAl29(t) J P 1 {r)P 2 {r + A/12) cxp 4 W w -^ (r+A/l2)) *(r)**(r + A/i 2 )dr + CC 



£ P JV(ii + t 2 ) + NVhhV(f 12 )S 12 ( f 12 ) exp 



-N^hT 2 V{-f 12 )S* 12 {-f 12 ) exp' 



-iA 12 e(t) 



(A.20) 



where A 12 9(t) = 9\{t) — 6 2 (t) is the phase delay between the two beams which samples the fringe pattern. In the Fourier 
plane, the interferometric equation takes the form: 



I{y) = Z P N(ti + t 2 ) + NVht 2 V(f 12 )S 12 (f 12 )S Ul2 + N^T 2 V*(f 12 )S* 12 (f 12 )5- 



(A.21) 



5^ is the Dirac function, and vi 2 is the temporal frequency at which the interferogram is modulated, that is the frequency 
of the moving piezoelectric mirror. 

Note that, at the difference of the image plane recombination, pupil plane interferometry transmits only the frequency 
baseline fi 2 . Nonetheless, the coherent flux, estimated in the Fourier plane at the frequency mirror (y = vi 2 ) takes the 
same form as in the image plane case, namely: 



F c 12 = NVhhV 12 S 12 (f 12 ) 



12 

and similarly the total number of detected photoevents is: 

K = I{0) = ZpN(h+t 2 ) 
and the normalized spectrum of the interferogram at the frequency ^12 is: 

/ ht 2 Vi 2 



i(y 12 ) 



(ti+t 2 )£. 



-Si 2 (f\ 2 ) 



(A.22) 



(A.23) 



(A.24) 



A. 3. 3. The interferometric phase error 

In previous sections, we have shown that independently of the chosen recombination scheme, the tot al flux, the cohcr - 
ent flux and consequently the normalized spectrum of t he interferogram can take the form of Eqs. (A. 14 A.16| A. 17), 
respectively. Using these expressions, Eqs. ( A.6| A.7 |A.8 ) can be rewritten as following: 



E(|Q| 2 ) = A^if 2 |U 12 | 2 < |Si 2 (/i2)| 2 >* +ZpN(h + t 2 ) + N pix a 2 det 
E(Q 2 ) = N 2 tit 2 V^ 2 < S 2 12 (f 12 ) 
E(Q) = NVhhV 12 < S 12 (f 12 ) > 



(A.25) 
(A.26) 
(A.27) 



NB: I Chelli (1989) has shown that the error 0} the phase does not depends on the object phase, hence one can consider in 
the following that the object is centro- symmetric, that is V\ 2 = \V\ 2 \ = Ke[Vi 2 }. 

It remains to derive in previous equations the moments of the interferometric transfer function S± 2 (fi 2 ): 



< S 12 (f 12 ) >* = J Pi(r)P2(r + A/ 12 ) < exp^^-^r+Xf^)) >$ dr 
= J P 1 (r)P2(r + A/i2)exp-2 < IW( r ^^( r+A/l2 »l 2>4 'dr 



jp exp 



(A.28) 
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< Sl 2 (fl 2 ) ><i 



Pi(r)P 2 {r + Xfi2)Pi{r')P 2 {r' + A/ 12 ) < eyi ^l(r)-^{r+*fi2)+^{r')~V 2 (v'+\f 12 )) >$ drdr ,/ 
Pi{r)P 2 {r + Xfi2)Pi{r')P 2 (r' + A/12) exp-*«** W-^(r-+A/ 12 )+^(r')-0 2 r (r'+A/ 12 ))| 2 >„ drd/ 
knowing that < 4 l i( r ) ( / ) i( r ') >$ = — ^^ r (r,r'), it comes, 

exp" 4 < y P 1 (r)P 2 (r + \fi 2 )Pi{r')P 2 {r' + \f 12 )eyi^' D *'^ ry ^' D ^ r+xfl2y+xfl2 ' ) drdr' 
changing axes reference: P%(r), P 2 (r) — > P(r) centered on 0, 



exp 



P(r)P(r')exp5 



P(r)P(r')exp3 



drdr' 



(A.29) 



< I^i2(/i 2 )| 2 >• = / ^lW^aCr + A/ 12 )P 1 (r')P 2 (r' + A/ 12 ) < e x P i («W-«(r+A/ la )-*r(r')+«(r'+A/ ia )) >$ drdr / 
= /" Pi(r)P 2 (r + A/i 2 )Pi(r')P 2 (r' + A/12) exp-5 < l^iW-^( r+A/l2 )-^( r ') + ^( r ' +A/ i2))l 2 >* drdr ' 
= y Pi(r)P 2 (r + A/i 2 )Pi(r')P 2 (r' + A/12) exp-^^VI^^+A/^y+A/^) drdr ' 



changing axes reference: 

P(r)P(r f )exp-* T '+<- ( - r ' r " , > 



P(r)P(r')exp-5 p *^ , '' r ') 



drdr' 



(A.30) 



Putting Eqs. (A. 28 A.29 A.30) into Eqs. (A. 25 A. 26 A. 27 1, we finally find that the variance of the phase can be written 
as the quadratic sum or 3 terms, corresponding respectively to the detector, photon and atmospheric regimes a"i = 



Vdet* + Vphot* + a a*m*> With: 



1 



det d 



2Y, P NH 1 t 2 \V 12 \ 2 e- 2a * 



l J \p(r)P(r')e-i' D i-r(ry)'] 2 drdr' - e" 4 ^, / \p{r)P{r')e^ v ^ r ' r ' 



drdr' 



[J P 2 (r)dr] e 



2 „-2o 



(A.31) 
(A.32) 

(A.33) 



We note that the expression of the variance in the photon noise case {pi hot ) is the generalization of the formula found 



by Vannier et al. (2006) (Eq. (14)), when the visibility is corrupted by the partially AO-corrected turbulent phase. 



A. 4. Phase noise in single-mode interferometry 

Because of the spatial filtering specific properties of single-mode devices, the single-mode interferometric equation changes 
slightly with respect to that of the multimode as we have to take into account the coupling efficiency which extenuates 
the number of coherent and incoherent photoevents. The coupling coefficients, respective ly for the photometric (pi(V), 
p 2 (V)) and interferometric (p 12 (V)) channels write ( |Mege et al.||2003| |Tatulli et al.|[2004| : 

n/ s (V*Si)(0) , h 9n 
Pi[V) = po r en, , . (i=[l,2J) 



Po 



P 12 (V) = p 



Po 



I V(f)Pi(r)Pi(r + A/) exp l «'^ T ')-^>+ A -> : » drd/ 
[JPi(r)dr] 2 
(V*S 12 )(f 12 ) 

J V{f - /i 2 )Pi(r)P 2 (r + A/) e ^Wi(r)-^(r+xf)) drd/ 
/ Pi(r)dr J P 2 (r)dr 



(A.34) 



(A.35) 
(A.36) 

where S®(f), S°,(f) are the photometric transfer functions of a perfect (a tmosphere-free) interferometer, and po the 
maximum coupling efficiency, shown to be ~ 80% ( Shaklan fc Roddierj|1988 ). 
And the single-mode interferometric equation writes: 

/(/) = VpNhpiQOHiif) + Z P Nt 2P2 (V)H 2 (f) + ^pNVhi2P 12 (V)H 12 (f) + zZ P NVW 2 (p 12 {V))* H* 12 {-f) (A.37) 
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where Hi(f), H 2 (f) and Hi 2 (f) are the normalized (i.e. -ffi(O) = PT 2 (0) = 1, and Hi 2 (fi 2 ) = 1) single-mode photometric 
and interferometric transfer functions, so-called the carrying wave (Mege 2002) which are fixed by the geometry of the 



single-mode device and therefore independent of the atmosphere. As a result, the single- mode coherent flux is defined 
by: 

F{ 2 = ^ P N^hhp 12 (V) (A.38) 

A. 4.1. Compact sources case 

In the case of compact sources - that is unresolved by the telescope, which is the most common case in interferometry - 
Equations |A.34| and A.35| are taking the simplified following form: 



, . f PAr)Pi{r + A/)exp l ^( r )-^ r ( r+A -f»drd/ , k . 

PitY) = Po J - ' { {[ P N - 2 = PoSi (A.39) 

[JPi(r)dr] 

12 J P 1 (r)P 2 (r + Xf) exp^M-Mr+V)) drdf 

p (V) = p V 12 - r r = PoV 12 pi2 A.40) 

J Pi(r)dr J P 2 (r)dr 

(A.41) 

where Si is by definition the instantaneous Strehl ratio relative to the i th telescope, and with: 

J P\{r)P 2 (r + A/) exp i ^i( r )~^( I '+ A /)) drdf 

Pl2 - /P 1 (r)dr/P 2 (r)dr (A ' 42) 

| pi2 1 2 being by analogy the instantaneous single-mode interferometric Strehl ratio. The interferometric equation finally 
rewrites: 

/(/) = HpNttpoStH^f) + X P Nt 2 p S 2 H 2 {f) + ^ P NVhf 2 V 12Pl2 H 12 {f) + ZpNVhhPoV^H^-f) (A.43) 
And the single-mode coherent flux, at the spatial frequency fi 2 takes the form: 

Fi 2 = Z P NVhhV 12 p 12 (A.44) 

whereas the number of photoevents writes: 

K = I(0) = H pNpoSi + £ P Nt 2Po S 2 ( A.45) 

which depends on the turbulent atmosphere through the fluctuations of the instantaneous Strehl ratio Si and S 2 . And 
the normalized spectrum of the interferogram is: 

i(/ia ) = VhhP^ (A.46) 

{J12> tiSi+t 2 S 2 v ' 

The generic equations of the phase noise thus rewrites: 

E(|Q| 2 ) = T,pN 2 tit 2 p 2 \Vi 2 \ 2 < \pi 2 \ 2 >* +Z P Nt lPo < Si >$ +Z P Nt 2 p < S 2 >$ +N pix a 2 det (A.47) 

E(Q 2 ) = E 2 P N 2 tit 2 p 2 V 2 2 < p\ 2 >* (A.48) 

E(Q) = ZpNVhhPoV^ < pi 2 >* (A.49) 
It now remains to compute the moments of the coupling coefficients: 

<S l >$ = S~ h (i=[l,2]) (A.50) 

where Si is the long exposure Strehl ratio of the i th telescope. Assuming equivalent AO correction for both telescopes, 
we have Si = S 2 = S 

J Pi(r)P 2 (r + Xf) < exp'WW-^fr+A/)) >$ drd/ _ j fi(r)^(r + A/)drd/ _ ff 2 

Pl2 * /Pi(r)dr/P 2 (r)dr 6 " / Pi(r)dr / P 2 (r)dr 6 " ( ' 

2 

which shows that the long exposure interferometric Strehl ratio is equal to the coherent energy e - °*>- . 

2 J Pi{r)P 2 (r + Xf)Pi(r')P 2 (r' + Xf) < expH««+«M-02>+A/)-02>'+A.f')) >$ drdr'd/d/' 

[J Pi(r)dr J P 2 (r)dr] 2 



-4a 



2 JPi(r)P 2 (r + Xf)Pi(r')P 2 (r' + Xf) exp5 p ^( r < r ')+5 x V(* , + A ./>'+ A /') drdr'd/d/' 



[J Pi(r)dr J P 2 {r)drY 
change of variables r + Xf — > u, r' + A/' — > v! : 
-4a 2 I Pi(r)Pi(r')exp^' D 't"-( r ' r "> drdr' / P 2 (r)P 2 (r') exps^^-'"') drdr' 

[/P 1 (r)dr/P 2 (r)dr] 2 



< \Pl2\ >* 
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JPi(r)P 2 (r + Xf)P 1 (r')P 2 (r' + Xf) < exp l W (r)-4>l(r')-V 2 (r+\f)+4> r 2 (r'+\f)) >$ drdr'dfdf 

[J P^dr J P 2 (r)dr] 2 
JP 1 (r)P 2 (r + Xf)P 1 (r')P 2 {r' + Xf) exp-'^ir^-i-D^r+Xfy+xf) drdr'd/d/' 

[f Pi(r)dr / P 2 {r)dr] 2 
change of variables r + Xf — > m, r' + Xf — > v! : 

J P 1 (r)P 1 (r')exp-^ v <t>^ r ' r "> drdr' f P 2 {r)P 2 (r') ey^~^ v ^ ry ^ drdr' 
[f Pi(r)dr / P 2 (r)dr] 2 



(A.53) 



Putting Eqs. ( |A.50||A.5l")|A.52[|Al)3"| into Eqs. ( |A.47||A.48}|A.49[ ), we find again that the variance of the phase can be 
written as the quadratic sum of 3 terms, corresponding respectively to the detector, photon and atmospheric regimes 



2 

"detj, 



2 -2< 



2 _ i (h + h)s 



2p E P M 1 t 2 |F 12 | 2 e- 2CT 5. 



atm^ 



f P(r)P{r')e-i v *rW>drdr' 

[/P(r)dr]V 2 



/ P(r)P(r')e3 x '^( r, >'-')drdr' 



(A.54) 
(A.55) 

(A.56) 



Appendix B: Coherent flux dropouts - Rician density probability 



By analogy with the study of Canales & Cagigal (1999) on the statistics of the instantaneous Strehl ratio in partial 
AO correction, the instantaneous interferometric Strehl ratio, namely / = |Ti 2 | 2 and I — \pi 2 \ 2 for multimode and 
single-mode interferometers respectively are shown to follow Rician density probability statistics of the form: 



dp(I) 



2a 2 



2a l J \ o~ 



(B.l) 



In this appendix, we aim to derive the theoretical exp ression of the moments of t he Rician distribution (namely a and a) 



in both multimode and single-mode cases. Following [Canales fc Cagigal (19991, the quantity a and a can be expressed 
from the moments of the real (A r ) and imaginary (Ai) parts of the complex instantaneous interferometric Strehl ratio 
(i.e. such as / = Af + A 2 ): 



a 4 = <A r >% +2 < A r >| (af - af) - (af - af ) 2 
2a 2 = af + af+ < A r >| -a 2 

where a 2 =< A 2 >* - < A r >i and a 2 =< A 2 >* - < A, >l 



(B.2) 
(B.3) 



B.l. Multimode case 

In multimode interferometry we have / = |Ti2| 2 , hence: 

/ Pi(r)P 2 (r + Xf) cos(^(r) - 0$(r + A/))dr 



A r = Re(T 12 ) = 
Ai = Im(T 12 ) 



/ Pi(r)P 2 (r + Xf) sin(#(r) - <%{r + Xf))dr 



(B.4) 
(B.5) 



Following the same approach as in appendix A. 3 we compute the first and second order moments of these quantities. 
/ Pi{r)P 2 (r + Xf) < cos(#(r) - $(r + Xf)) > $ dr 



< A r >$ = RefTi; 



exp <t>r 



„ , ^ t rrp \ [PlirWr + Af ) < sinffi(r-) - ^(r + Xf)) > $ dr 

< Ai >$ = Im(Ti 2 ) = - = 



(B.6) 
(B.7) 
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/ Pi(r)P 2 (r + Xf 12 )Pi(r')P 2 (r' + A/ 12 ) < cos(#(r) - <%(r + A/12)) coflfaftr') - $(r' + A/ 12 )) drdr' 



£ 2 



using cos(a) cos(&) = (cos(a — b) + cos(a + b))/2 

1 J Pi{r)P 2 {r + Xfx 2 )Px(r')P 2 (r' + A/12) < cos(#(r) - <^(r + A/12) - 0I(r') + 4> r 2 (r> + A/12)) >$ drdr' 



Y 2 



1 / Pi(r)P 2 (r + Xfi 2 )Pi(r')P 2 (r' + A/12) < cos(#(r) - 0J(r + A/12) + fl(r') - ^(r' + A/12)) >* drdr' 



changing axes reference: Pi, P 2 centered on 0, 

1 J [p(r)P(r')exp~5 I V( r < r ')| \p(r)P(r') exp-^r^': 



drdr' 



exp 



[fP 2 (r)dr] 2 

J [p(r)P(r')ex P 3 x '«r('':'-')l \p( r )P{r') exps 23 ^ ( r > r "> 
[JP 2 (r)dr] 2 



drdr' 



(B.8) 



< Aj >$ 



/ Pi(r)P 2 (r + A/ 12 )P 1 (r')P 2 (r' + A/ 12 ) < sin(#(r) - + A/ 12 )) sm(<^(r') 



^(r' + A/ia)) >* drdr' 



V 2 
2, p 



using sin(a) sin(6) = (cos (a — b) — cos (a + b) )/2 

1 J Pi(r)P 2 (r + A/ 12 )P!(r')P 2 (r' + A/ 12 ) < cos (^(r) - $(r + A/12) - <«(r') 
2 



^(r' + A/ 12 )) >$drdr' 



1 J P!(r)P 2 (r + A/i 2 )Pi(r')P 2 (r' + A/12) < cos (ff (r) - 4> r 2 {r + A/12) + #(r') 
2 

changing axes reference: Pi, P 2 centered on 0, 



^(r' + A/12)) >• drdr' 



X J P(r)P(T-')exp-2 z? ^<>> T1 ') 
2 



P(r)P(r')exp-2 c «,-('-^') 



drdr' 



[JP 2 (r)dr] 



-4<rl 



/ P(r)P(r')exp3 



P(r)P(r') exp 



drdr' 



[/P a (r)dr]' 



(B.9) 



B.6 


B.7 


B.8 


B.9 



in Eqs. (B.2 B.3|, we obtain the theoretical expressions of the parameters a and a of the 



timode CcLS6, ctS summarized in Table |B.1| of the appendix 



6.2. Single-mode case 

In single-mode interferometry we have / = |pi 2 | 2 , hence: 



A r = Re(p 12 ) = 
Ai = Im(pi 2 ) = 



$Pi(r)P 2 {r + A/) cos(#(r) - <%(r + A/))drd/ 

/ Pi(r)dr / P 2 (r)dr 
J Pi(r)P 2 (r + A/) Bin(0I(r) - #(r + A/))drd/ 



/ Pi(r)dr / P 2 (r)dr 

Following the same formalism as previously, we compute the first and second order moments of these quantities. 
/ Pi(r)P 2 (r + A/) < cos(#(r) - $(r + A/)) >4 , drd/ 



(B.10) 
(B.ll) 



< A r >$ 



/ Pi(r)dr / P 2 (r)dr 
/Pi(r)P 2 (r + A/) <sin(^(r)-^(r + A/)) > $ drd/ 
/Pi(r)dr / P 2 (r)dr 



= exp 
= 



(B.12) 
(B.13) 



< A 2 . >$ = 



< A 2 >$ 
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J P 1 (r)P 2 (r + A/)P(r')P 2 (r' + Xf) < cos(fl[(r) - <&(r + A/)) cos(^Qr') - ^ + A/')) > $ drdr'd/d/' 

[/P 1 (r)dT-/P 2 (T-)dr] 2 
using cos(a) cos(&) = (cos(a — b) + cos(a + 6))/2 

1 J Pi(r)P 2 (r + A/)P(r')P 2 (r' + Xf) < cosj^jr) - cb r 2 (r + A/) - #(r') + <%(r' + A/') >g drdr'd/d/' 

2 [/Pi(r)dr/P 2 (r)dr] 2 

| 1 J P(r)P 2 (r + Xf)P 1 (r')P 2 (r' + A/') < cos(^(r) - ^(r + A/) + #(rQ - <^(r' + A/')) > $ drdr'd/d/' 

2 [/P 1 (r)dr/P 2 (r)dr] 2 
change of variables r + A/ — > it, r' + Xf' — > u' : 

1 J P 1 (r)P 1 (r')exp-5 x '^('"' T '')drdr'/ P 2 {r)P 2 (r') exp-i^V^r-') drdr' 

2 [/P 1 (r)dr/P 2 (r)dr] 2 

^= f Pi (r)Pi(r')ex P 5 x '^('-^') drdr' f P 2 (r)P 2 (r') exps 27 ^,*-') drdr' 
+e 4<T *rJ ^ — i — i — — (B.14) 

[J P 1 (r)dr f P 2 (r)dr] 2 

J Pi(r)P 2 (r + Xf)P 1 (r')P 2 (r' + Xf) < smj^jr) - <j>\(v + Xf)) smj^jr') ~ + Xf')) > $ drdr'd/d/' 

[fp(r)dr fp>(r)dr] 2 
using sin(a) sin(6) = (cos(a — b) — cos(a + b))/2 

1 J Pi(r)P 2 (r + A/)P(r')P 2 (r' + A/') < cos(fl(r) - ^(r + A/) - 0£(rQ + <%(r' + A/') > $ drdr'df d/' 

2 [/P 1 (r)dr/P 2 (r)dr] 2 

1 f P(r)P 2 (r + A/)P(r')P 2 (r' + A/') < cos(0£(r) - ^(r + Af) + g(rQ - <^(r' + A/') > $ drdr'd/d/' 

2 [/P 1 (r)dr/P 2 (r)dr] 2 
change of variables r + A/ — >• u, r' + Xf' —> u' : 

1 J P^^P^rOexp-s^^^^drdr'/ P 2 (r)P 2 (r') exp-5 p *^ , '' r ') drdr' 

2 [/P 1 (r)dr/P 2 (r)dr 



2 



_ 4(T 2 / Pi(r)Pi(r') exps^Vr-KO drdr' / P 2 (r)P 2 (r') expi 2 ^ (*-,»•') drdr' 



, e -4^,. J + v / + v > r J ^ ' ^ ' r : (B.15) 

[J Pi(r)dr J P 2 (r)dr] 2 



B.13 B.14 B.15) in Eqs. ( |B.2 B.3), we obtain the theoretical expressions of the parameters a and a 



Putting Eqs. ( |B.12 

of the Rician distribution is the single-mode case, as summarized in Table |B.1| of the appendix. 

Appendix C: Astrometric phase error 

The astrometric phase error writes: 

al^H = 2*1 - 2cov(0 s , <t>? ah ) (C.l) 

where <f> s , and <f) Aah are respectively the partially AO-corrected phase of the astrophysical target and the phase estimated 
from the off-axis reference source. Act is the angular separation between both sources, and h is the altitude of the turbulent 
layer. 



Since ai has been computed in previous sections, it remains to compute the covariance part of Eq. (C.l). The formal 



expression of the covariance is given by Chelli (1989): 



_ 1 Re[E(QQ* Aah )] - Re[V{QQ Aah )} 
C ° V(0S '^ )_ 2 E(Q)E(Q Aq/i ) (C ' 2) 

where QaqH is the estimator of the coherent flux for the off-axis reference source. 
C.l. Multimode case 

Using formalism of Sect. [XJ one can derive the moments of the multimode estimator: 

E(Q) = N^hf 2 V 12 < S 12 (f 12 ) >4> (C.3) 

HQ Aah) = N Aah Vhhv£ ah < S Aah (f 12 ) > (C.4) 

HQQlah) = NN Aah txt 2 Vi2(V£ ah y < S 12 {f 12 ){S Aah )*(f 12 ) ><t> (C.5) 

E(QQAah) = NN Aah t 1 t 2 V 12 V 1 Aah < S 12 (f 12 )S Aah (f 12 ) >rP (C.6) 
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Table B.l. Parameters a and a of the Rician probability distribution function for both multimode and single- mode 
cases. 

Multimode case: 



exp 



2 — exp 



/ [p(r)P(r')es I V (T '' I '' ) ] 2 drdr' 



[JP 2 (r)dr] 



° = 2 



Single-mode case: 



1 [ S\P{r)P{r')c-^ T 't-^ r ' r ' ) ^ drdr' 



[/P 2 (r)dr] 



a 2 — exp 2<T * 



2 — exp 



[/P(r)P(r')es I, «- (7, ' r '' ) drdr'j' 



[/P(r)dr] 



° = 2 



1 [ [/ P(r)P(r')e-i v <t-r(ry) drAr >j' 



[JP(r)dr] 



where N Aah and V^ ah are respectively the number of photons (per surface unit and per time unit) and the visibility of 
the reference source, and with: 

< S Aah (f 12 ) >$ = J P 1 (r)P 2 (r + A/i 2 ) < exp i ^ ah W-^ ofc ( r+A/ia » >$ dr = £ P cxp-< (C.7) 

Since we assume one turbulent layer located at the altitude h, we have by definition (j)f ah (r) = </>i(r + Aah) and 
4>2 ah ( r ) = 4>2{r + Aah), hence: 



< S? 2 ah {f 12 ) >$ = J Pi(r)P 2 (r + A/12) < exp l ^ r+ * a V-^ r+Xfl2+ * ah » >$ dr = Epexp"^ 



(C.8) 



Furthermore, we have: 
< S\ 2 (f '\ 2 )Si 2 ah '(/ '12) >$ 



= J P 1 (r)P 2 (r + Xf 12 )Pi(r')P 2 {r' + A/ i2 ) < e xp J ^ 1 M-^( r+A/l2 )+^( r ' +Aa ' l )-^( T '' +A/l2+Aa ' l » > $ drdr' 
= J Pi(r)P 2 (r + Xf 12 )P 1 (r')P 2 (r' + A/12) e xp-5<l(^('-)-^^+^^)+^('-'+ Aa ' l )-^('-'+ A /^+ A «' l ))l 2 >* drdr' 
knowing that < <j>i(r)<j)i(r' + Aah) >$= cr^ — ^p^ r (r,r' + Aah), it comes, 

= cxp- 4 <- y Pi(r)P 2 (r + A/ 12 )Pi(r')P 2 (r' + A/ i2 ) exp kV4 r (ry+± a h)+iv^( r +\f 12 y+xf 12 +± a h) drdr , 
changing axes reference: Pi, P 2 centered on 0, 



exp 



P(r)P(r') Q-xphV+rirS+^hj] r P ( r )p( r ') cxp iX)^(rV+A«h) 



drdr' 



(C.9) 



< Si2(/i2)(S& a T(/i 2 ) >* (CIO) 
= y Pi(r)P 2 (r + A/i 2 )Pi(r / )P 2 (r / + A/12) < exp l ^ 1 M-^( r+A/l2 )-^( r ' +Aa/l ' + ^( T '' +A - fl2+Aa ' l » >$ drdr' 

= y Pi(r)P 2 (r + A/ 12 )P 1 (r')P 2 (r' + A/12) CX p-5 < l^ 1 ( T -)-^( r '+ A ^ 12 )-^( T -'+ Aa ' l )+^('-'+ A -' !l2 + Aa ' l ))l 2> * drdr' 
= y Pi(r)P 2 (r + A/i 2 )Pi(r')P 2 (r' + A/12) exp-s 1? +r('-.'''+Aa/.)-§i7* r (r+A/ 12 ,r'+A/ 12 +Aa/.) drdr ^ 
changing axes reference: Pi, P 2 centered on 0, 



y [p(r)P(r'))cxp-5 I, ^ (r ' r ' +Aa ' l) ] [p(rP(r') exp" 5 x, *- (r ' r ' +Aa ' l) 



drdr' 



(C.ll) 
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Putting Eqs. (C.7 C.8 C.9 C.lll into Eqs. (C.5 C.6 C.3 C.4), we finally find that the covariance part of the multimode 
astrometric phase writes: 



n 2 



cov(cj> s ,tf ah ) 



1 f P{r)P(r')e-^ v ^ {r ' r ' + ^ ahS > drdr' - e" 4CT «^ / P{r)P{r')e^ v *^ r ^ '+ Aah ) 



drdr' 



[/ P 2 (r)dr]V 2CT ? 



(C.12) 



C.2. Single-mode case 

We once again use formalism of Sect. [S] to compute the moments of the single-mode estimator: 
E(Q) = E P Ny/ht 2 ~p V 12 < P12 >4> 

E(Q Aah ) = z P N Aah Vhhp v£ ah < P ? 2 ah > <t> 

E(QQ Aah ) = Y ip NN Aah t 1 t 2P lV l2 v£ ah < p 12 pf 2 ah > 
with 

< P12 ><t> = e- a *r 
< p£ 2 ah >$ - c-< 
We also have: 



(C.13) 
(C.14) 
(C.15) 
(C.16) 

(C.17) 
(C.18) 



Aah Pi{r)P 2 (r + Xf)Pi(r')P 2 (r' + Xf) < c ^(Mr)+Mr'+^h)-<j, 2 (r+xf)-Mr'+xf+^ a h)) >$ drdr'd/d/' 

< P12P19 >$ = o 

[/Pi(r)dr/P 2 (r)dr] 2 

-Aal I Pi(r)P 2 (r + Xf)Pi(r')P 2 (r' + Xf) ^^^^y+^^+^v^ir+M-r'+M'+^och) drdr < d f d f 

[J Pi(r)dr J P 2 (r)dr] 2 
change of variables r + Xf — > u, r' + Xf' — > u' : 
_ 4(T 2 j P 1 {r)P 1 {r')e^ v ^ {ry+ * ah) drdr 1 / P 2 (r)P 2 (r') exp5 I5 ^ r > r '+ Aa ' 1 ) drdr' 

[/P 1 (r)dr/P 2 (r)dr] 2 ' '' J ' 

Pi(r)P 2 (r + Xf)P 1 {r')P 2 {r l + Xf) < e ^(Mr)-Mr'+^h)-Mr+^f)+Mr'+\f'+^h)) >$ drdr'd/d/' 

[J P 1 {r)dr J P 2 (r)dr] 2 

/ Pi(r)P 2 (r + Xf)Pi(r')P 2 (r' + Xf) e X p-| p ^(^^'+Ao./»)-|P^(r+Af,r-'+A/'+Ao ! fe) drdr ' d / d /' 

[/P 1 (r)dr/P 2 (r)dr] 2 
change of variables r + Xf —> u,r' + Xf' —> u' : 

j Pi (r)Px (r') exp- (r,r'+A*h) drdr > J p 2 ( r )p 2 (y ) exp - \T3^ T (r,r'+ Aa h) drdr , 



< pi2(p£r h r >* 



[J P 1 (r)dr f P 2 (r)dr] 2 

Putting Eqs. |C^[Cl8}[aT9}[a20] ) into Eqs. ( [ClIl[Cl^[ai3l[al4| , we find that the covariance part of the single-mode 
astrometric phase writes: 



(C.20) 



iAah\ 



J P{r)P(r')c-^ v *^ r ' r '+ Aah Mrdr'] - e~ i<r lr \j P(r)P(r , )ei v *r ( - r ' r ' +AoA '>drdr f 



[JP(r)dr]\ 



(C.21) 



Appendix D: Phase spatial fluctuations in partial AO correction 

The spatial fluctuations of the atmospheric phase can be characterized by the computation of the so-called spatial 
structure function 2?^(r,r') defined as the variance of the phase difference taken at two radii r and r': 



XV(r,r') = V^(r - r') =< \<f>(r) - <f>(r')\ 2 > 



(D.l) 



In the fully turbulent case , and assuming Kolmogorov's description of the atmosphere, the phase structure function can 
be written ( |Roddier||1979[ ): 

/ \ 

X>*(p) = 6.88(^J (D.2) 



E. Tatulli et al.: Multimode vs. single-mode interferometric phase, Online Material p 11 



tq being the fried parameter, that gives the ultimate spatial resolution achievable in presence of atmospheric turbulence, 
that is 9 tU rb = A/r . 

When partial AO correction is applied, low (spatial) orders of the corrugated wavefront are real-time compensated 
thanks to a deformable mirror. The variance cr? of the residual turbulent phase 4> r ( r ) not corrected by the system, can 



be approximated as following, providing a full correction of the first J< 1 Zernike modes of the turbulence Noll (1976): 

5/3 



0.2944J" 



-VV2 



Furthermore, the phase structure function of the residual phase can be written: 

< <jf{r) - 4> r (r + p) > 



1 



(D.3) 



(D.4) 



For p bigger than typically the distance between the mirror's actuators, the phases are uncorrelated, that is < <fi r (r) 
<j) r (r + p) >= 0, and consequently saturates at 2a\ . 

The long exposure AO corrected transfer function of the atmosphere is by definition: 



B(f) = cxp 



(D.5) 



Taking into account the properties of the phase structure function of the residual phase, it comes (Conan 1994) that 
£>(/) can be decomposed in two parts: one low frequency term - namely the halo {FTOhaio) ~ and one high frequency 

_ 2 

term saturating at e a <t>r , namely the coherent energy. It comes that the long exposure transfer function of the telescope 
T(f) — T (f).B(f), where T (/) is the transfer function of a perfect telescope, can be in first approximation modelled 
by a sum of two components: 



T(f) 



W).[e- 
To(f).[e- 



FTO halo (f)\ 
(1 - exp^*.-' 



exp 



(D.6) 



Eq. (D.4) also demonstrates that for interferometers equipped with AO systems, the phases 4>\ and 4> r 2 of telescopes 1 



and 2 are always uncorrelated. As a consequence, the phase structure function of the residual piston phase <^ 12 = <p x — <p 2 
is twice the phase structure function of the residual phase over one single telescope (assuming equivalent AO systems for 
both telescopes), that is T)^ 12 (p) = 2T> c f ir (p) . As a matter of fact, even if no AO correction is applied, the phases <j>i and 
4>2 can still be consider uncorrelated, as soon as the baseline is longer than the outer scale of the atmosphe re Cp. This is 
for ex ample the case for the VLTI at Paranal, where the outer scale has been estimated to be Co ~ 22m (Martin et al. 
20001. 



